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Magnetic properties of the natural mineral Cu2(OH)2C03 (malachite) were investigated through 
DC and AC susceptibility measurements. The analysis of the low-temperature part reveals a quan- 
tum spin-gap behavior with A ~ ISO-ftT. Consistently with the crystal structure, the magnetic 
susceptibility can be accurately described by a model of alternating chain. The non- frustrating 
residual inter-chain magnetic couplings, describing a sort of dimerized square planar lattice, are not 
strong enough to push the system towards a long-range ordered ground state, in good agreement 
with recent theoretical studies. 
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An intense renewal of interest for low dimensional 
quantum spin systems has stemmed from the Haldanc's 
conjecture [|l[, which pointed out the essential differ- 
ence between the ground state properties of the anti- 
ferromagnetic Heisenberg (AFH) chain with integer and 
half-integer spin. Quantum fluctuation may prevent low- 
dimensional AFH systems from falling in a long-range or- 
dered (LRO) Nel-like ground state for particular topolo- 
gies of magnetic interactions and spin values. For in- 
stance, Haldane proposed that the uniform spin 1 chain 
would exhibit a quantum disordered (spin liquid) ground 
state (Haldane phase) with a gap in the spin excitation 
spectrum, whereas the uniform spin 1/2 chain would be 
gapless and quasi-LRO. There has been considerable ev- 
idence in support of this statement during the last two 
decades on both experimental |0] and theoretical sides [pi . 
Interestingly, the uniform AFH spin chain is in a criti- 
cal state, in the sense that an inflnitesimal perturbation 
may push the system into a completely different ground 
state. In the case of exchange dimerization, the resulting 
system is an alternating chain which has a gapped dis- 
ordered ground state with a string order closely related 
with the Haldane phase's one jj]. 

One of the current issues in the field of quantum mag- 
netism is the occurrence of the spin liquid state in two- 
dimensional (2D) lattices, with or without the presence 
of magnetic frustration. In non-frustrated 2D AFH sys- 
tems, the case of dimerized square planar lattices, built 
up from alternating chain connected by inter-chain cou- 
pling, is especially interesting. With a vanishing dimer- 
ization within the chain and a strong enough inter-chain 
coupling, the uniform square planar, which ground state 
is LRO, can be recovered. On the other hand, a zero 
inter-chain coupling and a finite dimerization leads to a 
simple alternating chain, with a disordered spin liquid 
ground state. A quantum critical line therefore sepa- 
rates the two different ground states in the dimerization 
inter-chain coupling sub-space. Several theoretical stud- 
ies have been devoted to this subject 



In this context, we studied the magnetic properties of 
the malachite Cu2(OH)2C03. This well-known mineral, 
also used in jewelry, has, to our knowledge, never been 
investigated for its quantum magnetic properties. The 
crystal structure is monoclinic (space group P12i/al) 
with cefl parameters a=9.502 A, 6=11.974 A, c=3.240 
Aand 0=98.75° Q. The structure is built up from cor- 
rugated copper oxides (ac) planes formed with edge and 
corner-sharing distorted CuOg octahedra (see Fig. 1). 
The CO3 groups bridge the neighboring planes together 
and ensure the three-dimensional stability. The copper 
electronic configuration is 3cP (Cm^"*", s=1/2). In this 
paper, we show that the dimensionality of the magnetic 
interaction lattice defined by the different magnetic cou- 
plings between the spins 1/2 is intermediate between ID 
and 2D in Cu2(OH)2C03. Taking into account the four 
principal magnetic interactions, the structure appears to 
be quasi-2D consisting, within the ac plane, in a set of 
parallel alternating chains connected by weak inter-chain 
coupling, resulting in a kind of dimerized square planar 
lattice. Despite the inter-chain couplings, the magnetic 
susceptibility can be accurately described using the high 
precision fitting procedure recently developed by John- 
ston and co-workers |10|| for an quantum disordered al- 
ternating chain. This clearly indicates that inter-chain 
couplings are too weak to push the system towards a 
LRO ground state, as expected from recent theoretical 
studies 



I. EXPERIMENTAL 

Two different natural samples of Cu2(OH)2C03, orig- 
inating one from the Urals, Russia and one from Zaire, 
were used in this study. Their purity was checked by 
X-rays diffraction and Electron Probe Microscope Anal- 
ysis. The samples were found to be single phased and no 
impurities were detected. DC and AC magnetic measure- 
ments were performed in a commercial Quantum Design 
MPMS-5S SQUID magnetometer. The fields were set to 



1 kOe for the DC measurement and to 2 Oe for the AC 
one, with a driving frequency of 997 Hz. 

The rather complex fitting formula for the magnetic 
susceptibility of an alternating chain proposed by John- 
ston and co-workers [lOJ was inserted in a program in- 
cluding a Simplex fitting algorithm. We checked the 
absence of any error among the eighty four constants 
required in the formula by fitting successfully numeri- 
cal data obtained for alternating chains by the transfer- 
matrix density-matrix renormalization group (TMRG) 
method llOJ for a =0 (dimer), 0.5, 0.95 and 1 (uniform 
spin 1/2 chain). 



II. RESULTS 

Figure 2 shows the magnetic susceptibility data mea- 
sured at 1 kOe. The most striking feature is the presence 
of a broad maximum in the xC^) curve at 120 K. This 
behavior is typical either of low dimensional antiferro- 
magnets or of systems with a gap in the spin excitation 
spectrum. The drop of the susceptibility below the broad 
maximum is balanced below 25 K by a Curie-like tail. 
In order to detect a possible spin gap, we performed a 
low temperature fit which describes the susceptibility as 
the sum of a T-independent Xo term plus a Curie- Weiss 
C/{T — 9) law and a thermally activated susceptibility 
^7^-i/2g-A/T g^itable for most of the quasi-lD spin gap 



systems 11 1. The good quality of the fit, depicted in the 
inset of Fig. 2, strongly supports the existence of a spin 
gap in the malachite. The low temperature tail corre- 
sponds to 2.4 percent of free spin per Cu. The value 
of the spin gap loosely depends on the fitted tempera- 
ture range (typically 2-60 K) and stay within the 110- 
130 K range. Interestingly, it roughly corresponds to the 
position of the maximum of the susceptibility T{xmax)^ 
implying a ratio A/T(xmax) close to unity. As a con- 
sequence, the malachite is probably not a trivial collec- 
tion of isolated dimers of spin 1/2, for which the ratio 
A/T{xmax) is around 1.6. We also attempted to fit our 
data using the expression for isolated dimers, with the 
number of dimers as a free parameter. For any reason- 
able value of the g factor (2 < g < 2.25), it leads to an 
unrealistic value of number of dimers, between 60 and 70 
percent of the expected value. 

As we will see below, magnetic frustration could pos- 
sibly exist in the malachite. We therefore measured the 
AC susceptibility Xac=x'+^x" of the malachite to detect 
a possible spin glass-like behavior at low temperature. It 
results (i) in a perfect matching between the component 
x' and the DC susceptibility and (ii) in the absence of 
any peak in the x" part down to 2 K, consistent with the 
absence of any ordered state down to low temperature. 



III. DISCUSSION 

In order to get more insight into the real magnetic 
interaction network of the malachite, it is important to 
establish at least a hierarchy of the different magnetic 
couplings existing in this system. The leading magnetic 
coupling here originates in Cu-O-Cu superexchange. The 
main parameters affecting the strength and the sign of 
superexchange are the Cu-0 distance, the Cu-O-Cu an- 
gle and the crystal structure via the Madelung poten- 
tial |12|. The Anderson-Kanamori-Goodenough rules [l^ ] 
state that the super-exchange Cu-O-Cu is the sum of a 
weak angle-independent ferromagnetic coupling and of a 
strong antiferromagnetic coupling maximum for (p— 180° 
and tending to zero at 90°. The total magnetic inter- 
action is therefore expected to cross zero and to change 
its sign close to 90°. A theoretical approach based on a 
three bands Hubbard model and dedicated to the study 
of several cuprates with edge-sharing CUO4 units showed 
that this crossing occurs close to = 95° |1J] . 

The above description of the structure in terms of cor- 
rugated planes of octahedra is somewhat misleading since 
most of these octahedra have a broad distribution of 
copper-oxygen distances ranging from 1.89 Ato 2.67 A 
Q. By restricting the pertinent Cu-0 distances to the 
1.89-2.11 Arange, all the polyhedra change from CuOg 
octahedron to distorted CUO4 squares planar and the 
2D connectivity is lost within the ac plane, as shown in 
Fig. 3(a). The effective magnetic network only consists 
in parallel chains running in the ac plane and by only 
two types of magnetic interactions Ji and J2 alternating 
along the chains (see Fig. 3(b)). A detailed description 
of Ji and J2 is given in Table I. By comparison with the 
phase BaCu2Si207 [^, where almost similar Cu-O-Cu 
bonds lead to J=280 K, a rough estimate gives 200 < Ji 
< 400 K. As the couphng J2 is mainly mediated through 
a Cu-O-Cu bond with an angle of 106° (122° for Ji), h 
must be significantly smaller than Ji. 

By decreasing order of intensity, the next magnetic 
couplings should be located within the ac planes and link 
the chains together. Three new magnetic inter-chain cou- 
plings J3, J4 and J5 (see Fig. 3(c) and Table I) can be dis- 
tinguished and are all mediated via Cu-O distances close 
to 2.37 A. J5 should be negligible compared to J3 and 
J4, since the involved Cu-O-Cu angles, 97.5° and 97.7°, 
are very close to the critical value (95° in the cuprates) 
for which the magnetic superexchange vanishes |lj] . The 
network including the four main interactions Ji, J2, J3 
and J4, shown in Fig. 3(c), has strong similarities with 
the non-frustrated dimerized square lattice treated theo- 
retically by Singh and co-workers [5] and depicted in Fig. 
3(e). Estimating the J3 and J4 values is not a simple 
task; the Cu-O-Cu angles in J3 and J4 are however close 
to the one appearing in Ji and J2, respectively (see Table 
I). Therefore the main difference originates in the Cu-0 



distances (2.37 Aagainst 1.91 A). Assuming a power-law 
dependence r~" of the superexchange strength with dis- 
tance, a value n=10-ll can be inferred from the work of 
Mizuno et al. |l|] on Cu-0 networks. As (2.37/1.91)"" 
~ 0.1, J3 and J4 can be estimated to be roughly one order 
of magnitude smaller that Ji and J2, respectively. Note 
that, contrary to J3 and J4, the interaction J5 could lead 
to a magnetic frustration since it couples spins over the 
diagonal of the dimerized square. 

The last kind of magnetic coupling that can be consid- 
ered are those mediated by the CO3 groups. It results in 
(i) a next nearest neighbor (NNN) interaction within the 
J1-J2 chain and (ii) in a set of inter-plane couplings J/p, 
inducing a strong frustration as shown in Fig. 3(d). Fi- 
nally, Dzyaloshinsky-Moriya interactions could possibly 
exist in the malachite because of the absence of inversion 
center between the different neighboring Cu ||l^. They 
are however most probably negligible since the magni- 
tude of the low temperature susceptibility is hardly rec- 
oncilable with the existence of a weak ferromagnetism. 

From the above magneto-structural analysis, it ap- 
pears that a simple model of alternating chains is a good 
starting point to understand the magnetism of the mala- 
chite. As previously discussed, the ratio /S./T{xmax) is 
close to one in the malachite. In an alternating chain 
where strong bonds J alternate with weaker bonds aJ 
(ajl), such a feature is expected for a close to 0.5 Q. 
Also, as T{xmax)/ J loosely depends on the dimcrization 
a in an alternating chain {T{xmax)/J ^ 0.6 Ijl^Dj the de- 
duced value of the strongest magnetic coupling within the 
chain in the malachite is Ji ~ 200 K. This corresponds to 
the lower limit of Ji estimated from magneto-structural 
arguments. To be more quantitative, we tried to fit the 
susceptibility curve of our two samples using a fitting 
procedure for alternating chains recently proposed in ref. 
[10] and suitable for any dimerization parameter < a 
< 1. We used the formula: 



X{T) = Xo + C/iT - 



+Ng^fI%/{kBJ) * Xalt.chain{T, J, Of) 



(1) 



where Xo is a T-independent term, C the Curie con- 
stant, 9 the Weiss temperature and J the strongest cou- 
pling within the alternating chain. The free parameter 
were Xo, C, 9, g, J and a. The best result of the fit, 
depicted in Fig. 2, was obtained for g = 2.08 (2.20), J = 
196K (203 K) and a = 0.49 (0.53). The values in brackets 
correspond to the results obtained on the second sample 
(not shown in Fig. 2). The most important parameters 
of the fit, J and a, appear to be very similar from one 
sample to another. The resulting spin gap can be com- 
puted using the formula A/ J « {l+af/^ {l-af/'^ 0] 
and is equal to 131 (128) K. This value compares remark- 
ably well with the one obtained from the low-temperature 
approach (A=126 K). 

The good agreement between data and the global fit 



proposed by Johnston and co-workers supports the model 
of simple alternating chains. However, a fit of the mag- 
netic susceptibility is not sufficient to prove that the mag- 
netic network considered is correct. It is therefore im- 
portant to check if a quantum spin-gap state is expected 
taking into account the presence of (i) the possible NNN 
coupling within the chain, (ii) the intcr-chain coupling J3 
and J4. 

The coupling induced by the CO3 groups can safely 
assumed to be small. It would otherwise result in a 3D 
highly frustrated magnetic lattice, which is hardly recon- 
cilable with the existence of a large spin gap in the mala- 
chite. As far as the inter-chain couplings J3 and J4 are 
concerned, it is possible to be more quantitative. Singh 
and co-workers [51 considered the topology depicted in 
Fig. 3(e), similar as the dimerized square planar de- 
scription of the malachite but with J2=J3=J4=J'. They 
found a quantum critical point for ac=(J'/J)c=0.39 ± 
0.01, separating the LRO ground state (a > ac) from 
the spin liquid disordered one (a < ac)- Qualitatively, it 
is expected that ac will tend to the critical value of a pure 
alternating chain (ac=l) in the limit of small inter-chain 
(J3,J4) coupling, which is the case in the malachite. This 
intuition is confirmed by the work of Katoh and Imada 
[Q • They considered a slightly different dimerized square 
planar topology shown in Fig. 3(f), where strong bonds 
J alternate with weak bonds aJ . As the inter-chain cou- 
pling J3=J4=J_L decreases from J_l/J=0.70 to 0.26, the 
critical dimerization ac increases from 0.40 to 0.74. 

In the malachite, our previous estimate of the ratio 
intcr-chain over intra-chain couplings, Jj^/J « J3/J1, and 
of the dimerization parameter a were J^/J « 0.1 and a 
K, 0.50, respectively. From the above discussion, it can 
be inferred that the critical dimerization Uc for the mala- 
chite is larger than a « 0.50 for an inter-chain coupling 
J^/J « 0.1. A disordered spin liquid ground state is 
therefore expected in the malachite, in agreement with 
the low temperature fit of the susceptibility. 



IV. CONCLUSION 

In summary, the natural mineral Cu2(OII)2C03 ex- 
hibits a quantum spin liquid ground state, consistent 
with its structure including s=l/2 alternating chain with 
rather low inter-chain couplings. Also, the malachite pro- 
vides a good example of interesting magnetic properties 
in a natural mineral. In the field of quantum magnetism, 
other copper-oxide based minerals could possibly be a 
rich source of new materials that has only been poorly 
investigated so far. 
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FIG. 3. (a) Copper-oxygen lattice obtained by keeping only 
the Cu-0 distances below 2.11 A.(b) Equivalent magnetic lat- 
tice, where magnetic couplings are represented as solid and 
dashed segments, (c) Magnetic lattice including the four lead- 
ing couplings, (d) Details of the interactions within (along the 
a axis) and perpendicular to (along the b axis) the chains, in- 
cluding those possibly mediated by the COa groups (see text), 
(e) and (f ) Magnetic topologies considered in Ref. [5] and Ref. 
[7], respectively. 



FIG. 1. Crystal structure of the malachite Cu2(OH)2C03. 
The main structural units are the CuOe octahedra and the 
CO3 triangles. 



FIG. 2. DC magnetic susceptibility of Cu2(OH)2C03. 
Data are depicted in open squares, the fit of the suscepti- 
bility following the model of alternating chains as a solid line 
and the intrinsic susceptibility Xmeasured — Xo — C jiT — 9) 
as a dashed line. The inset shows the same data in a dou- 
ble-logarithm plot as well as the low temperature fit, as de- 
scribed in the text. 




Figure 1, E. Janod et al. 
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Figure 2, E. Janod et al 
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Figure 3, E. Janod et al. 



Table 



Table I : classification by decreasing order of strength of the different magnetic interactions in 
the malachite and their description including the Cu-0 distances and Cu-O-Cu angles. Several 
cuprates systems from the literature are also included for comparison. A "number of Cu-O-Cu 
paths involved" of one (two) indicates comer (edge) -sharing CUO4 or CuOg units. 



Number of 
Cu-O-Cu 

paths 
involved 



dcui-o do- 



Cu2 



Cu-O-Cu angle 



Estimate of J's 



1 



1.92 A 1.90 A 



122° 



1.91 A 
2.05 A 



1.91 A 
2.11 A 



106° 
94.7° 



1.90 A 2.37 A 



118° 



1.91 A 2.37 A 



106° 



1.91 A 

1.92 A 



2.37 A 
2.37 A 



97.7° 
97.5° 



BaCuzSizOv [16] 


2 


1.96 A 


1.97 A 


124° 


280 K 


BaCuzGezO? [16] 


2 


- 


- 


135° 


540 K 


CuGe03 [15] 


2 


1.93 A 


1.93 A 


99° 


150 K 



